ON A DOUBLY NONLINEAR DIFFUSION MODEL OF CHEMOTAXIS 
WITH PREVENTION OF OVERCROWDING 
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QQ Abstract. This paper addresses the existence and regularity of weak solutions for a fully parabolic model 

of chemotaxis, with prevention of overcrowding, that degenerates in a two-sided fashion, including an extra 
nonlinearity represented by a p-Laplacian diffusion term. To prove the existence of weak solutions, a Schauder 
fixed-point argument is applied to a regularized problem and the compactness method is used to pass to the 
limit. The local Holder regularity of weak solutions is established using the method of intrinsic scaling. The 
results are a contribution to showing, qualitatively, to what extent the properties of the classical Keller-Segel 
chemotaxis models are preserved in a more general setting. Some numerical examples illustrate the model. 

< 

,— , 1. Introduction 

1.1. Scope. It is the purpose of this paper to study the existence and regularity of weak solutions of 
the following parabolic system, which is a generalization of the well-known Keller-Segel model [ , 2, 3] of 
chemotaxis: 

d t u-div(\VA(u)\ p - 2 VA(u)) + div ( X uf(u)Vv) = in Q T := SI x (0, T), T > 0, SI C R N , (1.1a) 
i ^ i dtv — dAv = g(u, v) in Qt, (1.1b) 

\VA(u)\ p - 2 a(u)—- = 0, — =0 on S T :=aOx(0,T), (1.1c) 
orj orj 

u(x, 0) = Uo(x), v(x,0) = vo(x) on (l.ld) 

On 

where Q C R is a bounded domain with a sufficiently smooth boundary dfl and outer unit normal rj. 
Equation (1.1a) is doubly nonlinear, since we apply the p-Laplacian diffusion operator, where we assume 
2 < p < oo, to the integrated diffusion function A(u) := J Q n a(s) ds, where a(-) is a non-negative integrable 
function with support on the interval [0, 1]. 

In the biological phenomenon described by (1.1), the quantity u = u(x,t) is the density of organisms, 
such as bacteria or cells. The conservation PDE (1.1a) incorporates two competing mechanisms, namely 
the density-dependent diffusive motion of the cells, described by the doubly nonlinear diffusion term, and a 
motion in response to and towards the gradient \7v of the concentration v = v(x,t) of a substance called 
chemoattractant. The movement in response to Vv also involves the density-dependent probability f(u(x,t)) 
for a cell located at (x, t) to find space in a neighboring location, and a constant \ describing chemotactic 
sensitivity. On the other hand, the PDE (1.1b) describes the diffusion of the chemoattractant, where d > 
is a diffusion constant and the function g(u,v) describes the rates of production and degradation of the 
chemoattractant; we here adopt the common choice 

g{u,v) = au — f3v, a,f3>0. (1-2) 
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We assume that there exists a maximal population density of cells u m such that f(u m ) = 0. This 
corresponds to a switch to repulsion at high densities, known as prevention of overcrowding, volume-filling 
effect or density control (see [ ]). It means that cells stop to accumulate at a given point of Q after their 
density attains a certain threshold value, and the chemotactic cross-diffusion term \ u f{ u ) vanishes identically 
when u > u m . We also assume that the diffusion coefficient a(u) vanishes at and ^ m , so that (1.1a) 
degenerates for u = and u = u m , while a(u) > for < u < u m . A typical example is a(u) = eu(l — u m ), 
e > 0. Normalizing variables by u = u/u m , v = v and f{u) = f(uu m ), we have u m = 1; in the sequel we will 
omit tildes in the notation. 

The main intention of the present work is to address the question of the regularity of weak solutions, which 
is a delicate analytical issue since the structure of equation (1.1a) combines a degeneracy of p-Laplacian type 
with a two-sided point degeneracy in the diffusive term. We prove the local Holder continuity of the weak 
solutions of (1.1) using the method of intrinsic scaling (see [ , ]). The novelty lies in tackling the two types 
of degeneracy simultaneously and finding the right geometric setting for the concrete structure of the PDE. 
The resulting analysis combines the technique used by Urbano [ ] to study the case of a diffusion coefficient 
a(u) that decays like a power at both degeneracy points (with p = 2) with the technique by Porzio and 
Vespri [ ] to study the p-Laplacian, with a(u) degenerating at only one side. We recover both results as 
particular cases of the one studied here. To our knowledge, the p-Laplacian is a new ingredient in chemotaxis 
models, so we also include a few numerical examples that illustrate the behavior of solutions of (1.1) for 
p > 2, compared with solutions to the standard case p = 2, but including nonlinear diffusion. 

1.2. Related work. To put this paper in the proper perspective, we recall that the Keller-Segel model is a 
widely studied topic, see e.g. Murray [ ] for a general background and Horstmann [ ] for a fairly complete 
survey on the Keller-Segel model and the variants that have been proposed. Nonlinear diffusion equations for 
biological populations that degenerate at least for u = were proposed in the 1970s by Gurney and Nisbet 
[ ] and Gurtin and McCamy [10]; more recent works include those by Witelski [ ], Dkhil [ ], Burger et al. 
[13] and Bendahmane et al. [ ]. Furthermore, well-posedness results for these kinds of models include, for 
example, the existence of radial solutions exhibiting chemotactic collapse [14], the local-in-time existence, 
uniqueness and positivity of classical solutions, and results on their blow-up behavior [ ], and existence 
and uniqueness using the abstract theory developed in [16], see [17]. Burger et al. [ ] prove the global 
existence and uniqueness of the Cauchy problem in M N for linear and nonlinear diffusion with prevention of 
overcrowding. The model proposed herein exhibits an even higher degree of nonlinearity, and offers further 
possibilities to describe chemotactic movement; for example, one could imagine that the cells or bacteria are 
actually placed in a medium with a non-Newtonian rheology. In fact, the evolution p-Laplacian equation 
u t = div (| Vu\ p ~ 2 Vu), p > 1, is also called non-Newtonian filtration equation, see [ ] and [ , Chapter 2] 
for surveys. Coming back to the Keller-Segel model, we also mention that another effort to endow this 
model with a more general diffusion mechanism has recently been made by Biler and Wu [ ], who consider 
fractional diffusion. 

Various results on the Holder regularity of weak solutions to quasilinear parabolic systems are based 
on the work of DiBenedetto [ ]; the present article also contributes to this direction. Specifically for a 
chemotaxis model, Bendahmane, Karlsen, and Urbano [4] proved the existence and Holder regularity of 
weak solutions for a version of (1.1) for p = 2. For a detailed description of the intrinsic scaling method and 
some applications we refer to the books [5, 6]. 

Concerning uniqueness of solution, the presence of a nonlinear degenerate diffusion term and a nonlinear 
transport term represents a disadvantage and we could not obtain the uniqueness of a weak solution. This 
contrasts with the results by Burger et al. [13], where the authors prove uniqueness of solutions for a 
degenerate parabolic-elliptic system set in an unbounded domain, using a method which relies on a continuous 
dependence estimate from [21], that does not apply to our problem because it is difficult to bound Av in 
L°°(Qt) clue to the parabolic nature of (1.1b). 

1.3. Weak solutions and statement of main results. Before stating our main results, we give the 
definition of a weak solution to (1.1), and recall the notion of certain functional spaces. We denote by p' 
the conjugate exponent of p (we will restrict ourselves to the degenerate case p > 2): - + \ = 1. Moreover, 
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C w (0,T,L 2 (tt)) denotes the space of continuous functions with values in (a closed ball of) L 2 (Q) endowed 
with the weak topology, and (•, •) is the duality pairing between W 1}P (Q) and its dual (W 1,P (Q)Y . 

Definition 1.1. A weak solution of (1.1) is a pair (u,v) of functions satisfying the following conditions: 

< u(x,i) < 1 and v(x,i) > for a.e. (x,i) G Qt, 

ueC w (0, T, L 2 (^)) , dtu G L p/ (0, T; (W 1 *^))') , u(0) = u , 

A(u) = [ a(s) ds G L p (0, T; W 1 ' P (Q)) , 
Jo 

v G L°°(Q T ) H U (0, T; jy x ' r (SI)) H C(0, T, L r (^)) for all r > 1, 

and, for all <p G L^(0, T; W 1 '^)) and ^ G L 2 (0, T; H 1 ^)), 

(dtu,<p)dt + J J ^\VA(u)\ p - 2 VA(u)- X uf{u)VvYv^dxdt = {) 1 

I (dtv , ip) dt + d II Vv-Vipdxdt= II g(u^v)ip dx dt. 

Jo JJqt JJqt 

To ensure, in particular, that all terms and coefficients are sufficiently smooth for this definition to make 
sense, we require that / G C^O, 1] and /(l) = 0, and assume that the diffusion coefficient a(-) has the 
following properties: a G C^O, 1], a(0) = a(l) = 0, and a(s) > for < s < 1. Moreover, we assume that 
there exist constants S G (0, 1/2) and 72 > 71 > 1 such that 

71000 < a ( s ) < 72^(5) for sG [0,(5], 71^(1 - s) < a(s) < 72 ^(1 - 5) for s G [1 — J, 1] , (1.3) 
where we define the functions 0(s) := s^/b-i) and ^(s) := s^/b-i) f or (3 2 > f3\> 0. 
Our first main result is the following existence theorem for weak solutions. 

Theorem 1.1. If u$,vq G L°°(Q) with < uo < 1 and > a.e. in Q, then there exists a weak solution 
to the degenerate system (1.1) in the sense of Definition 1.1. 

In Section 2, we first prove the existence of solutions to a regularized version of (1.1) by applying the 
Schauder fixed-point theorem. The regularization basically consists in replacing the degenerate diffusion 
coefficient a{u) by the regularized, strictly positive diffusion coefficient a £ (u) := a(u) +e, where e > is the 
regularization parameter. Once the regularized problem is solved, we send the regularization parameter e to 
zero to produce a weak solution of the original system (1.1) as the limit of a sequence of such approximate 
solutions. Convergence is proved by means of a priori estimates and compactness arguments. 

We denote by dtQr the parabolic boundary of Qt, define M := |M|oo,q t , and recall the definition of the 
intrinsic parabolic p-distance from a compact set K C Qt to dtQr as 

p-dist(K; d t Q T ) := inf (\x-y\ + M {p ~ 2)/p \t - s\ 1/p ). 

(x,t)eK, (y,s)ed t QT 

Our second main result is the interior local Holder regularity of weak solutions. 

Theorem 1.2. Let u be a bounded local weak solution of (1.1) in the sense of Definition 1.1, and M = 
Halloo, q t - Then u is locally Holder continuous in Qt, i.e., there exist constants 7 > 1 and a G (0,1), 
depending only on the data, such that, for every compact K C Qt, 

1, \ / \ 1 ^A\xi-x 2 \+M( p - 2 V p \t 2 -t 1 \ 1 { p Y w/ w N Tjr 

{ p-dist(K-dtQT) J 

In Section 3, we prove Theorem 1.2 using the method of intrinsic scaling. This technique is based on 
analyzing the underlying PDE in a geometry dictated by its own degenerate structure, that amounts, roughly 
speaking, to accommodate its degeneracies. This is achieved by rescaling the standard parabolic cylinders 
by a factor that depends on the particular form of the degeneracies and on the oscillation of the solution, and 
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which allows for a recovery of homogeneity. The crucial point is the proper choice of the intrinsic geometry 
which, in the case studied here, needs to take into account the p-Laplacian structure of the diffusion term, 
as well as the fact that the diffusion coefficient a(u) vanishes at u = and u = 1. At the core of the proof is 
the study of an alternative, now a standard type of argument [ ]. In either case the conclusion is that when 
going from a rescaled cylinder into a smaller one, the oscillation of the solution decreases in a way that can 
be quantified. 

In the statement of Theorem 1.2 and its proof, we focus on the interior regularity of w, that of v follows 
from classical theory of parabolic PDEs [ ]. Moreover, standard adaptations of the method are sufficient 
to extend the results to the parabolic boundary, see [5, 23]. 

1.4. Outline. The remainder of the paper is organized as follows: Section 2 deals with the general proof of 
our first main result (Theorem 1.1). Section 2.1 is devoted to the detailed proof of existence of solutions to 
a non-degenerate problem; in Section 2.2 we state and prove a fixed-point-type lemma, and the conclusion 
of the proof of Theorem 1.1 is contained in Section 2.3. In Section 3 we use the method of intrinsic scaling 
to prove Theorem 1.2, establishing the Holder continuity of weak solutions to (1.1). Finally, in Section 4 
we present two numerical examples showing the effects of prevention of overcrowding and of including the 
p-Laplacian term, and in the Appendix we give further details about the numerical method used to treat 
the examples. 

2. Existence of solutions 

We first prove the existence of solutions to a non-degenerate, regularized version of problem (1.1), using 
the Schauder fixed-point theorem, and our approach closely follows that of [ ] . We define the following closed 
subset of the Banach space L p (Qt)- 

K := {u G L P (Q T ) : < u(x,i) < 1 for a.e. (x,i) G Qt}- 

2.1. Weak solution to a non-degenerate problem. We define the new diffusion term A £ (s) := A(s)+es, 
with a £ (s) = a(s) + e, and consider, for each fixed e > 0, the non-degenerate problem 

d t u e - div (\X7A £ (u £ )\ p - 2 VA £ (u £ )) + div ( X f(u £ )Wv £ ) = in Q T , (2.1a) 

d t v £ - dAv £ = g(u £1 v £ ) in Q T , (2.1b) 

\VA £ (u £ )r 2 a £ (u £ )^ =0, ^ = on E T , (2.1c) 

u £ {x, 0) = uo(x), v e (x, 0) = vo(x) for x G ft. (2. Id) 
With u G /C fixed, let v £ be the unique solution of the problem 

d t v £ — dAv £ = g(u,v £ ) in Q Tl (2.2a) 
dv 

= on E T , v e (x, 0) = v (x) for x e ft. (2.2b) 

Given the function v e , let u £ be the unique solution of the following quasilinear parabolic problem: 

d t u £ - div (\VA £ (u £ )\ p - 2 VA £ (u £ )) + div ( X u e f(u e )Vv e ) = in Q T , (2.3a) 

\VA £ (u £ )\ p ~ 2 a £ (u £ )^ = on E T , u £ (x,0) = u (x) for x G ft. (2.3b) 

Here vq and uq are functions satisfying the assumptions of Theorem 1.1. 

Since for any fixed u G /C, (2.2a) is uniformly parabolic, standard theory for parabolic equations [ ] 
immediately leads to the following lemma. 

Lemma 2.1. Ifv G L°°(Q), then problem (2.2) has a unique weak solution v £ G L°°(Q T )nL r (0, T; W 2 ' r (ft))fl 
C(0, T; L r (ft)) ; for all r > 1, satisfying in particular 

KIIl~(q t ) + \\ v s\\l^(o,T;L^(Q)) < C, \\v £ \\ L 2^ T . H i m < C, \\d t v E \\ L *(Q T ) < C, (2.4) 
where C > is a constant that depends only on \\vq\\ L oo^, ol, (3, and meas(QT)- 
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The following lemma (see [ ]) holds for the quasilinear problem (2.3). 

Lemma 2.2. If uq G L°°(^) ; then, for any e > 0, there exists a unique weak solution u £ G L°°(Qt) H 
LP(0 : T;W^ p (n)) to problem (2.3). 

2.2. The fixed-point method. We define a map O : /C — » /C such that 0(iZ) = where u £ solves (2.3), 
i.e., 6 is the solution operator of (2.3) associated with the coefficient u and the solution v £ coming from 
(2.2). By using the Schauder fixed-point theorem, we now prove that O has a fixed point. First, we need 
to show that © is continuous. Let {^ n }neN be a sequence in /C and u G K be such that u n — > u in L p (Qt) 
as n — > oo. Define i/ £n := 0(tZ n ), i.e., ii £n is the solution of (2.3) associated with iZ n and the solution v £n of 
(2.2). To show that u £n — > 0(u) in L p (Qt) 7 we start with the following lemma. 

Lemma 2.3. TTie solutions u £n to problem (2.3) satisfy 

(i) < u £n (x,t) < 1 /or a.e. (x, £) G Qt- 

(ii) The sequence {u £n } neN is bounded in L p (0, T; W 1,P (Q)) D L°°(0, T; L 2 (^)). 
(hi) TTie sequence {u £n } ne ^ is relatively compact in L P (Q T ). 

Proof The proof follows from that of Lemma 2.3 in [] if we take into account that {d t u £n } ne ^ is uniformly 
bounded in i/ (0,T; (W^ P (Q)Y). ^ □ 

The following lemma contains a classical result (see [22]). 

Lemma 2.4. There exists a function v £ G L 2 (0, T; i^ 1 (^)) such that the sequence {v £n } ne ^ converges 
strongly to v in L 2 (0, T; H 1 (ft)). 

Lemmas 2.2-2.4 imply that there exist ^ £ G L p (0,T; W 1,p (ft)) and v £ G L 2 (0, T; i^ x (ft)) such that, up to 
extracting subsequences if necessary, u £n — > ii £ strongly in L p (Qt) and v £n — > v £ strongly in L 2 (0,T; i7 1 (ft)) 
as n — > oo, so 6 is indeed continuous on JC. Moreover, due to Lemma 2.3, 0(/C) is bounded in the set 

W:={ueL p (0,T;W^ p (n)) : c^ii G L p ' (0, T; (W 1 ' p (ft)) / ) }. 

Similarly to the results of [ ], it can be shown that W L p (Qt) is compact, and thus is compact. Now, 
by the Schauder fixed point theorem, the operator has a fixed point u £ such that Q(u £ ) = u £ . This implies 
that there exists a solution (u e ,v £ ) of 

jT (d t u e ,<p) dt + JJ [\VA £ (u £ )\ p - 2 VA £ (u £ )- X u £ f(u £ )Vv £ y 

I (d t v £ ,ijj) dt + d I 1 Vv £ -V^dxdt= 1 1 g(u £ ,v £ )^ dx dt, (2.5) 

Jo JJqt JJqt 

V<p G L p (0,T; W hp {n)) and G L 2 (0,T; H 1 (Q)). 

2.3. Existence of weak solutions. We now pass to the limit e — » in solutions (ii £ ,v £ ) to obtain weak 
solutions of the original system (1.1). From the previous lemmas and considering (2.1b), we obtain the 
following result. 

Lemma 2.5. For each fixed e > 0, the weak solution (u £: v £ ) to (2.1) satisfies the maximum principle 

< u e (x,t) < 1 and v e (x,t) > for a.e. (x,t) G Q T - (2.6) 
Moreover, the first two estimates of (2.4) in Lemma 2.1 are independent of e. 

Lemma 2.5 implies that there exists a constant C > 0, which does not depend on e, such that 

\\ V £\\l°°(Q t ) + II^IIl 00 ^, T;L 2 (^)) - ^, \\V £ \\ L 2 (^T;// 1 (Q)) - ^ ' 

Notice that, from (2.6) and (2.7), the term g(u £ ,v £ ) is bounded. Thus, in light of classical results on L r 
regularity, there exists another constant C > 0, which is independent of £, such that 

\\ d t v e\\ L r(Q T ) + \\ V e\\ L r(0,T;W^ m < C for all T > 1. 
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Taking ip = A £ (u £ ) as a test function in (2.5) yields 

f (dtu e ,A(u e ))dt + e [ (d t u £ ,u £ )dt + ff \VA £ (u £ )\ p dx dt - ff X f(u £ )Vv £ -VA £ (u £ )dxdt = 0; 
Jo Jo JJqt JJqt 

then, using (2.7), the uniform L°° bound on u £l an application of Young's inequality to treat the term 
Vv £ ■ V ' A £ (u £ ), and defining A £ (s) := J Q S A £ (r) dr, we obtain 

sup / A £ (u £ )(x,t)dx + e sup / l^i^D dx+ ff |V 'A £ (u £ )\ p dx dt < C (2.8) 
o<t<rJn o<t<rJn 2 JJq t 

for some constant C > independent of e. 

Let (f G £ p (0, T; W llP (Q)). Using the weak formulation (2.5), (2.7) and (2.8), we may follow the reasoning 
in [4] to deduce the bound 

\\dtU£\\ LP ^o,T;(w^p(n)y) - (2.9) 
Therefore, from (2.7)-(2.9) and standard compactness results (see [24]), we can extract subsequences, which 
we do not relabel, such that, as e — > 0, 

f A e (u £ ) — > A(u) strongly in L p (Qt) and a.e., 
u £ — > u strongly in L q (QT) for all > 1, 
v £ — > v strongly in L 2 (Qt), 

Vv £ weakly in L 2 (Q T ) and VA £ (^ £ ) VA(w) weakly in L P (Q T ), (2.10) 

|VA e (u e )|*- 2 VA e (u e ) ri weakly in L p> (Q T ), 
v £ ^v weakly in L 2 (0, T; i/ 1 (Q)), 

,9 t u e <9 £ ?i weakly in i/ (0,T; (W 1 ^))') and <9 £ v weakly in L 2 (0,T; (i/ 1 ^))'). 

To establish the second convergence in (2.10), we have applied the dominated convergence theorem to 
u £ = A~ 1 (A £ (u £ )) (recall that A is monotone) and the weak-* convergence of u £ to u in L°°(Qt)- We also 
have the following lemma, see [4] for its proof. 

Lemma 2.6. The functions v £ converge strongly to v in L 2 (0, T; H 1 ^)) as e — > 0. 

Next, we identify Ti as |VA(ii)| p_2 V-A(i£) when passing to the limit e — > in (2.5). Due to this particular 
nonlinearity, we cannot employ the monotonicity argument used in [4]; rather, we will utilize a Minty-type 
argument [ ] and make repeated use of the following "weak chain rule" (see e.g. [ ] for a proof). 

Lemma 2.7. Let b : R — > R 6e Lipschitz continuous and nondecreasing. Assume u G L°°(Qt) is such that 
d t u G L^(0,T;(W 1 ^(ft)) / ), G Z7(0, T; W 1 ^)), u(a;,0) = u (x) a.e. on ft, ^ G L°°(fi). //we 
define B(u) = J^b(^)d^, then 

— / (d t u,b(u)(j)) dt = / / B(u)d t (j) dx dt + / B(uo)(j)(x,0) dx — / B(u(x, s))(/)(x, s) dx 
Jo Jo Jq Jq Jq 

holds for all G £>([0,T] x ft) and for any 5 G (0,T). 

Lemma 2.8. TAere AoZd I\ = | VA(^)|^- 2 VA(^) and V ' A £ (u £ ) -> VA(ix) strongly in L P (Q T ). 
Proof. We define Qt •= {(t, s,x) : (x, s) G Q^, t G [0,T]}. The first step will be to show that 

(Ti - \Va\ p ~ 2 Va) • (VA(u) - Va) dx ds dt > 0, V<j G 1^(0, T; W 1 '^)). (2.11) 

For all fixed £ > 0, we have the decomposition 

(\V A £ (u £ )\ p - 2 V A £ (u £ ) - \Vo-\ p - 2 V(j) ■ (VA(u) - Va) dxdsdt = h + / 2 +/ 3 , 

/i := /// |VA £ (^)|^ 2 VA £ (^) • (VA(tx) - VA £ (n £ )) dxdsdt, 
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h:= /// (\VA e {u e )\ p - 2 VA e (u e ) - |V<r| p - 2 V<r) • (VA £ (u £ ) - V<r) dxdsdt, 

J J J Q T 

I 3 := HI \Vcr\ p - 2 Vcr • (VA £ (u £ ) - VA(u)) dxdsdt. 



T 



Clearly, I2 > and from (2.10) we deduce that I3 — > as £ — > 0. For ii, if we multiply (2.1a) by 
G L p (0,T; W 1,:p (^)) and integrate over Qt, we obtain 



/ / (d t u £ ,(f))dsdt- HI xUef(u £ )Vv £ -V(j)dxdsdt+ 1 1 1 \V A £ {u £ )\ p ~ 2 V A £ {u £ ) -V (j) dx ds dt = {). 

JO JO J J J Q T J J J Qt 

Now, if we take <j> = A(u) - A £ (u £ ) G Z7(0,T; W llP (Q)) and use Lemma 2.7, we obtain 

h=- I I (d t u e ,A(u)) dsdt+ f / (0 t iz e , A e (u e )) dsdt 
JO t/ Jo Jo 

+ XUef(u £ )Vv £ - (VA(u) - VA £ (u £ )) dxdsdt 

J J J Q T 

= — / / (dtu e ,A(u))dsdt+ // A £ (u £ )dxdt — T / A £ {uo)dx 
Jo Jo JJqt J n 

+ /// XUsf(u £ )Vv £ • (VA{u) - VA £ {u £ )) dxdsdt. 

J J J Q T 

Therefore, using (2.10) and Lemma 2.6 and defining A(u) := J Q n A(s) ds, we conclude that 

lim 1\ = — / / (dtu, A(u)) ds dt + / / A(u(x , t)) dx dt — T / A(uo(x))dx, 
Jo Jo j^ Jq 

and from Lemma 2.7, this yields I\ — >• as £ — > 0. Consequently, we have shown that 

lim JJJ {\V A £ (u £ )\ p ~ 2 V A £ (u £ ) - \\/a\ p ~ 2 \/a) • (V A(u) - Va) dx ds dt > 0, 

which proves (2.11). Choosing a = A(u) - A£ with A G R and £ G L p (0,T; W liP (Q)) and combining the two 
inequalities arising from A > and A < 0, we obtain the first assertion of the lemma. The second assertion 
directly follows from (2.11). □ 

With the above convergences we are now able to pass to the limit e — > 0, and we can identify the limit 
(u,v) as a (weak) solution of (1.1). In fact, if ip G L p (0, T; W 1,P (Q)) is a test function for (2.5), then by 
(2.10) it is now clear that 

/ (d t u £ ,(p)dt^ / (d t u,(p) dt as e — > 0, 
Jo Jo 

II \VA £ (u £ )\ p - 2 VA £ (u £ ) -Vipdxdt -> II \VA(u)\ p - 2 VA(u)-Vcpdxdt as e -> 0. 
JJq t JJqt 

Since ft(u e ) = u £ f(u £ ) is bounded in L°°(Qt) and by Lemma 2.6, v £ — > u in L 2 (0,T; ff 1 ^)), it follows that 
// x u ef(u £ )X7v £ 'V(pdxdt^ // x u f( u )^ v '^ ( fdxdt as £ — > 0. 

JJq t JJqt 

We have thus identified u as the first component of a solution of (1.1). Using a similar argument, we can 
identify v as the second component of a solution. 
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3. Holder continuity of weak solutions 

3.1. Preliminaries. We start by recasting Definition 1.1 in a form that involves the Steklov average, defined 
for a function w G L 1 (Q T ) and < h < T by 

nt+h 



( i r + 

I - / w(-,r)dr if t G (0,T-ft], 
[O if t G (T- h,T}. 



Kx 



Definition 3.1. A local weak solution for (1.1) is a measurable function u such that, for every compact 
K C Q and for all < t < T - h, 

[d t (u h )ip + (\VA(u)\*- 2 VA(u)) h ■ V<p - ( X uf{u)Vv) h • V(^} dx = 0, Vy? G (3.1) 

{^} 

The following technical lemma on the geometric convergence of sequences (see e.g., [ , Lemma 4.2, Ch. I]) 
will be used later. 

Lemma 3.1. Let {X n } and {Z n }, n G No ; be sequences of positive real numbers satisfying 

x n+1 < Cb n (xi+ a + X^Z 1 ^) , Z n+1 < Cb n (x n + zi+ K ) , 

where C>1, b>l, a>0 and k > are given constants. Then X n , Z n — > as n ^ oo provided that 

X + < (2C)- (1+ * )Aj 6-( 1+/ ^ 2 , with a = min{a,«}. 

3.2. The rescaled cylinders. Let B p (xo) denote the ball of radius p centered at xq. Then, for a point 
(xo,to) ^ ^ n+1 , we denote the cylinder of radius p and height r by 

(xo,t ) + Q(t, p) := B p (x ) x (t - r, t ). 

Intrinsic scaling is based on measuring the oscillation of a solution in a family of nested and shrinking 
cylinders whose dimensions are related to the degeneracy of the underlying PDE. To implement this, we fix 
(xo,to) G Qt] after a translation, we may assume that (xo,£o) = (0>0)- Then let e > and let R > be 
small enough so that Q(R P ~ £ , 2R) C Qt, and define 

/i + := ess sup u, /i := essinf u, uo := essosc u = /i + — pT . 

Q(Rp-e,2R) Q(RP~z,2R) Q(RP—,2R) 

Now construct the cylinder Q(aoR p , i?), where 

'o;\2-p 1 

with m to be chosen later. To ensure that Q(aoR p , R) C Q(R P ~ £ , 2i?), we assume that 

*(*o <"» 

and therefore the relation 

essosc u <uo (3.3) 

Q(a RP,R) 

holds. Otherwise, the result is trivial as the oscillation is comparable to the radius. We mention that for uj 
small and for m > 1, the cylinder Q(aoR p , R) is long enough in the t— direction, so that we can accommodate 
the degeneracies of the problem. Without loss of generality, we will assume uo < S < 1/2. 
Consider now, inside Q(aoR p , R), smaller subcylinders of the form 

Q% = (0,t*) + Q(dRV,R), d =Q 2 ~ P m J 4)]P -i > **<°- 

These are contained in Q(a R p ,R) if a R p > —t* + dRP, which holds whenever 0(w/2 m ) < tj)(u>/A) and 

„ f(uj/2) 2 -PRP _ (uj/2)p- 2 Rp > 
1 € V iPiu/^P- 1 ~ <P(Lo/2 m )P- 1 ' , 



a 
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These particular definitions of ao and of d turn out to be the natural extensions to the case p > 2 of their 
counterparts in [ ]. Notice that for p = 2 and a(u) = 1, we recover the standard parabolic cylinders. 

The structure of the proof will be based on the analysis of the following alternative: either there is a 
cylinder Q f R where u is essentially away from its infimum, or such a cylinder can not be found and thus u is 
essentially away from its supremum in all cylinders of that type. Both cases lead to the conclusion that the 
essential oscillation of u within a smaller cylinder decreases by a factor that can be quantified, and which 
does not depend on uo. 

Remark 3.1. (See [8, Remark A.2)) Let us introduce quantities of the type BiR e uo~ hi , where Bi and bi > 
are constants that can be determined a priori from the data, independently of uo and R, and depending only 
on N and p. We assume without loss of generality, that 

BiR e cj- b * < 1. 

If this was not valid, then we would have uo < CR £ for the choices C = max^ B\ jh ande = 9/ min 7 - bi , and 
the result would be trivial. 

3.3. The first alternative. 

Lemma 3.2. There exists v$ G (0, 1) ; independent of uo and R, such that if 

|{0M)eQ£ : u(x,t) > 1 -w/2}\ < Vo\Qr\ (3.4) 
for some cylinder of the type Q f R , then u(x,t) < 1 — uo/A a.e. in Q^/ 2 * 
Proof. Let u u := min{^, 1 — cj/4}, take the cylinder for which (3.4) holds, define 

R R 

~2 + 2 n + i: 

and construct the family 



Rn 



n G N , 



QRn 



(0,**) + Q(dR p n1 R n ) = B Rn x (r n ,t*), r n := t* - dR p n1 n G N ; 



note that Q 1 



Rn 



!R/2 



oo. Let {£ n }neN be a sequence of piecewise smooth cutoff functions satisfying 



fn = 1 in 



2?i+i 2^( n + 1 ) 



( n = on the parabolic boundary of Q l R 

2p(n+l) 



(3.5) 



R 



dRP 



and define 



UO 

4 



n G N . 



Rp 



2n+2 ' 

Now take ip = [(u^h ~ ^n] + £n> K = Br u in (3.1) and integrate in time over (r n ,t) for t G (r n ,t*). Applying 
integration by parts to the first term gives 



F 1 := / / d 8 u h [(u u ) h - k n \ + ^dxds 

Jr n J B Rn 

= \f [ d s ({[(u u ) h -k n ] + ) 2 )&dxds+(l-^- 
= \\ (K - k n ]t) 2 e n dxds -\f (K - k n ]+) 2 e n dxds 

I I {K-k^+f^ds^ndxds 



1 + 



^ dx ds 



p 

2 
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In light of standard convergence properties of the Steklov average, we obtain 



*i - F* 



1 



B Rn x{t} 



V 



r n J B R 



( 1 -i- fc ")(l,x(j"-( 1 -i)] +S<iXdS 



p 



^ds^ndxds) as/i^O. 



lB Rn x{r n } 

Using (3.5) and the nonnegativity of the third term, we arrive at 



Rp 



X{uuj>k n } 



dx ds 



UJ ~i*~J T J BRn 



X{u>i-lj/4} dx ds 



T n JB Rn 



X{ Uuj >k n } dx ds, 



, B Rn x{t} 

the last inequality coming from u> 1 — uj/A > k n . Since [u u — k n ) + < cj/4, we know that 

(k, - k n } + ) 2 = (K - k n ] + f- p (K - k n ] + ) P > Q 2 ~ P (K - k n ] + Y > Q 2 ~ P (K - k n ] + ) p ; 

therefore, the definition of d implies that 



BR n x{t} 



T n J B Rn 



X{ Uuj >k n } dx ds. 



(3.6) 



We now deal with the diffusive term. The term 



converges for h — > to 



where we define 



F 2 := f I {a(uf- l \Vu\ p - 2 Vu) h ■ V{[( Uu ,) h - k n }+&} dxds 

Jr n J B Rn 

O 

K ■= f I aiuf-^Vur^u • (VK, - k n )+a+p(uu, ~ fc„)+Cr 1 V^) dxds 

Jr n J B Rn 

= f f a(u) p - 1 \£ n V(u u; -k n ) + \ p dxds + F*, 

Jr n J B Rn 

F* :=p [ [ a^uf-^Vu^Vu-V^u^ - fc*)^" 1 dxds. 



Since Viu^ — k n ) + is nonzero only within the set {k n < u < 1 — cj/4} and 

a{u) > 7i'0(cc;/4) on {k n < u < 1 — cj/4}, 
we may estimate the first term of F 2 * from below by 



r n JB Rn 



«H P_1 knV(^ - k n )+\ P dxds > [jMuj/4)} 



p-1 



^V(w w - &n) + \ P dx ds. (3.7) 



T n JB Rn 



Let us now focus on F 2 *. Using that V(w a; — k n ) + is nonzero only within the set {k n < u < 1 — a;/4}, 
integrating by parts, and using (1.3) and (3.5), we obtain 



K^rMvK - k n )+\ p 'iv^ik - fc„)+er i dxds 



Tn JB R „ 
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p [i---k r 



t„ JB Rn 



\ p-1 
a(s) ds) J> dx ds 

o/4 J + 



< P [72^(^/2)] 



p-1 



/ / iV^nK^o, - /c n ) + |^nV(^ a; - k n ) + \ P 1 dxds 

Jr n J B Rn 



+ P (j)~lj B {L\ /4 a(S)dS ) P + ((P-^n'^nf + en-^Qdxds 



Next, we take into account that 



and apply Young's inequality 



1-00/4 



a(s) ds I < —if) (cj/4) . 



ab<— a r + ifa,6>0, - + — = 1, e > 0, 



(3.8) 



for the choices 



r=p, a = |V£„|(w w - fc„) + , 6 = |V(u w - k n ) + \ p and e x p 



y ( 7 r - imov^ 

P 7 2 P "V(^/2)f- 1 



> 0. 



This leads to 



^2* < 



p _i /cj\p 2^ n+1 ) y* 



dx ds 



i^^r 1 (^) 

+ (p- l)e?' [72^/2)]^ T / |^nVK - kn) + \ p dxds 



-hp (^-J il){u/A) p Rp J X{u„>k n } dx ds 



< 



( 7 r X - 1)^(^/4) 



p _i (uy 2^ n+1 ) 



X{ Uuj >k n ] 



dx ds 



+ ( 7 r 1 -l)^(^/4) p - 1 T / \i n V{u„-k n )+\ p d X ds 

Jr n JB Rn 



2 / iJ \ p -1 2 p ( n+1 ) Z* 1 Z 4 

-hp (^-J il){u>/l) p Rp J X{ Uuj >k n }dxds. 



Hence, from (3.7) and (3.9), and observing that 



L^(w/4)J 



2 J 



we obtain 

F2>ti>(u>/*y 



I / |£nV(ua, - k n ) + \ P dxds 

Jr n J B Rn 



p 2 + 2 Pfr 



Finally, for the lower order term 



r vii 1 




[p(7r 1 -i)J 





-) Rp ^ (^/ 4 ) p y y^ x{u u >k n }dxds. 



F 3 := T / (x^/WV^ • V{[K)^ - £ n ] + ^}dx^ 

Jr n J B Rn 



(3.9) 



(3.10) 
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we have 



F 3 -> F 3 * := T / Vv • (V(« w - + - fc„)+er 1 Ve„) d^ds 



,B R n 

aXuf(u)Vv • V(^a; — & n ) + £n ^ ^ s 



> / / x w /( w ) • V^ n (w w - fc n ) + (n 1 dxds as ft, — > 0. 
Applying Young's inequality (3.8), with 



r = p, a = VK-g + $ n , ^X^/Her'Vv and eg = ^(o;/^" 1 > 0, 



using the fact that (u^ — k n ) + < uo/A and defining M := \\x u f( u )\\L°°(Q T )i we ma Y estimate F| as follows: 

P ^ JB Rn p'e P 2 J Tn JB Rn 

+ pM J r j B \ Vv \{^)\ V ^\X{u^>k n }dxds 

< Uiu/Af- 1 f [ \V(u^-k n )+^dxds+ {p/2) ~ P/P M ] f [ \Vv\v' X{Uuj>kn} dxds 
2 JtJb r J V' ^{u>/4)J Tri J B 



fuo\P / pM y l I ' 

7 / / \y£,n\ P X{ Uuj >k n }dxds + — -r / / | Vv\ p X{ Uuj >k n } dx ds, 

V4/ Jr n JB Rn p'tl Jr n JB Rn 



applying again Young's inequality (3.8) to the last term in the right-hand side, this time with 

r=p, a= |V^|c^/4, & = M|Vv|, eg = ^ (w/4) > 0. 

Using (3.5), we obtain 

F* <F** :=^(^/4f _1 T / \V{u u -k n ) + ^ n \ v dxds 

1 Jr n J B R 



MP 



-P'/P 



/ / \Wv\ p 'x {Uu) > kn} dxds 



+ (4) Rp i> (^/ 4 ) p y y^ X{« w >fc n }^d5- 

Additionally, using Holder's inequality, we may write 

f £ /* / / f f \ x ~ x Ip 

Jr J Br ^ P ' X ^^ dXdS ^^<^(Q T )[J T \ A LRM d ^) ' 

where ^ (<j)| denotes the measure of the set 

A k n ,R n ( a ) := { x e B Rn : u(x,a) > k n }. 

Thus we obtain 



Ft < 2^/4) P_1 y y^ |£„VK-fc„) + | P ^ s + (-) — — 1 p( j/4) p - 1 J X {u u >k n }dxds 

Q~ P ' /P +p] m\U lQT) (/ K,*>)l ^ 



p'tp (w/4) 



(3.11) 
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Combining the resulting estimates (3.6), (3.10), (3.11) and multiplying by 2(u/2) p 2 yields 
ess sup / ([uu — kn]^) P £n d% ds + - j j \^ n V (u u — k n )~^ \ P dx ds 

r n <t<t* JBtuxU} d Jr n Jb* 



< \ ^p2 p ~ 2 + p 2 + 2** 



P'l2 



Phr 1 - 1) 



P- 1 } L >/'<": " 2 



(u/2) p ~ 2 MP' 
p'tjj (cj/4) 



+ 2 



4 7 Bp d 

•t* 



X{u w >fe n } 



<ix (is 



(3.12) 



llv^llU 



Lp'p(Qt) 



1-1/p 



Next we perform a change in the time variable putting t := \{t — t*), which transforms Q(dR^, R n ) into 
Furthermore, if we define u UJ {-^ : E) := u^-^t) and £ n (*>£) = £n then defining for each n, 



A n := / Xiu^kr,} dx dt 

J-R p n JB Rn 

we may rewrite (3.12) more concisely as 



| (So, - fcn) + e"n||^ (Q ^ ) < 2 -^" 2 +p 2 + 2^ 



X{u w >fc n } ^ 



+ 2 



fp\-p'/v 

(f) +p 



P ( 7 f- 1 - 1) 

M p ' /w\(p-2)/p 



}(;)' 



(3.13) 



where V p (ft r ) = L°°(0, T; L p (f2)) fl L p (0, T; W 1,p (Cl)) endowed with the obvious norm. Next, observe that 
by application of a well-known embedding theorem (cf. [5, §1.3]), we get 



2p( 



n+1) V 



( ) A n+1 -\k n /c n +l| P ^n+l < ||(^w kn ^\\l,Q(R p n+1 ,R n+1 ) ~\\( U " kn ^ + ^ n \\ P p,Q{R p n ,R n 



4/ 



< elite - fc i+f ir * a p/{n+p) 



(3.14) 



Now, applying (3.13), we get 
1 foo\P A _ f 3 



2p(n+l) 



Q P A n+1 <2C^p2^+p" + 2^ 



+ 2C 



pItT 1 -!). 

M p ' /w\(p- 2 )/p 



) p2p( : +i) jl i + ,>/(N +> „ 



/p\-q/p ltf 1 /w\ 
(2) +^(2) 



^( W /4) 1 - p - 1 /P||Vt;|r L ' !) , p(QT) Ai- 1 /p+p/(^+p). 

(3.15) 



Now let us define 



Xr, 



\Q(Rn, R n ) 



A 1/p 



, n e N Q . 



Dividing (3.15) by ^ru (lY \Q(K+i,Rn+i)\ yields 



X n+1 < 2 pn ^2c| ^2 p " 2 + p 2 + 2 p02 



PI2 



(I) 



-p'/p 



LpW _1 -i)J 



j^l+p/^+p) 



< 7 2 p "(Xi + «+X:zi +K ), neN , 
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with a p/(N + p) > 0, ft = p - 2 > and 



7 := 2Cmax<J -p2 p ~ 2 + p 2 + 2 P/32 



/ v 

Pl2 



hhr 1 - 1) 



(i) 



-p7p 



^{2) [W4)] 



p-i 



v p-2 



> 0. 



(In the choice of ft we need the assumption that p is strictly larger than 2.) In the spirit of Remark 3.1, let 
us assume that 

Q P ~ 2 [^(co/4)r^R^<l. 

Therefore, with this assumption we conclude that 7 is independent of uo and R. 
Reasoning analogously, we obtain 

Z n+ i < j2 pn (X n + Z\ +K ) . 
Now, let a = min{a, ft} and notice that, if we set v$ := 27"( 1+/ ^ /cr (2 p )"( 1+/ ^ /cr2 , it follows from (3.4) that 

-(i+/^)/cr^2P)-(i+^)/c r2 > (3.16) 



X + Z x + K < 2 7~ 



Then, using Lemma 3.1, we are able to conclude that X n , Z n — > as n — > 00. Finally, notice that — > i?/2 
and fc n — > 1 — cj/4, and this implies that 

|{(M) e Q((R/2) p ,R/2) : ^(x,F) > 1-cj/4}| = \ {(x,t) G Q% 2 : ufot) > 1 - w/4} | = 0. 
This completes the proof. □ 

Now we show that the conclusion of Lemma 3.2 is valid in a full cylinder of the type Q(r,p). To this 
end, we exploit the fact that at the time level — i := t* — d(R/2) p , the function x 1— > u(x,i) is strictly below 
1 — uj/4 in the ball B R j 2 - We use this time level as an initial condition to make the conclusion of the lemma 
hold up to t = 0, eventually shrinking the ball. This requires the use of logarithmic estimates. 

Given constants a, 6, c with < c < a, we define the nonnegative function 



In- 



' a + c — (s — b)\± 
whose first derivative is given by 



In 



a + c ± (6 — s) 
if - < 



if 6±c^ s ^ 6± (a + c), 



5 ^ b ± c, 



(3.17) 



if 6±c^s^6± (a + c) 



< 



[0 ifs^6±c 
and its second derivative, away from s = b ± c, is 

fe )C )" = {(^Ac)'} 2 >0. 
Given it bounded in (xo,to) + Q(r, p) and a number fc, define 

H^k'= ess sup |(w-/c) ± |, 

(a5 ,to)+Q(r,p) 

and the function 

:= 4± fcifciC («), 0<c< J ff„ ± fe . 

Lemma 3.3. For even/ number v\ G (0, 1), £/iere exists s± G N, independent of uo and R, such that 
\{x G : u(x,i) > 1-uj/2 Si } \ < V!\B R/2 \ for all t G (-£,0). 



(3.18) 
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Proof. Let k = 1 — uo/A and 



V2 2+n , 



(3.19) 



with n G N to be chosen. Let < Q{x) < 1 be a piecewise smooth cutoff function defined on B R / 2 such that 
£ = 1 in B R /4 and |V£| < C/R. Now consider the weak formulation (3.1) with tp = 2g+ (uh)(g + )' '(uh)( p 
for K = B R j 2l where £> + is the function defined in (3.17). After an integration in time over (— t, t), with 
t G (— t, 0), we obtain Gi + G 2 — G3 = 0, where we define 

Gi := 2 / / d s {u h }g + (u h )(g + Y(u h )C p dxds, 

J —t J B R / 2 



G 2 
G 3 



t JB 



t JB 



(\VA(u)r 2 a(u)Vu) h ■ V{e + (u h )(Q+)'(u h )C p } dxds, 
( X uf(u)Vv) h ■ V{g + (u h )(g + Y(u h )C P } dxds. 



R/2 



Using the properties of the function £, we arrive at 
rt 



Gi= / / d s {g+(u h )} 2 edxds= f {Q+(u h )} 2 Cdx- f {g+ (u h )} 2 ( p dx. 

J-tJB R/2 JB R/2 x{t} JB R/2 x{-t} 

Due to Lemma 3.2, at time — t, the function x 1— > u(x,t) is strictly below 1 — uo/A in the ball B R / 2 , and 
therefore g+(u(x, —i)) = for x G Consequently, 

6?i-> / {g + (u)} 2 C p dx- I {q + (u)} 2 ( p dx = f {g+ (u)} 2 ( p dx as h 0. 

JB R/2 x{t} JB R/2 x{-t} JB R/2 x{t} 

(3.20) 

The definition of implies that 

u — k < H+ k ess sup ^ - 1 + ^ 



Q(t,R/2) 



(3.21) 



If = 0, the result is trivial; so we assume H^ k > and choose n large enough so that 



< 



22+n ^ JJ n,/c* 



Therefore, since -f/c — iz + c>0, the function g + {u) is defined in the whole cylinder Q(i, R/2) by 



In 



Kk+ c + k ~ u 



if > /c + c, 



otherwise. 



Relation (3.21) implies that 



ff+ fc +c+fc-u 2c 4 



< 4 - = 2 n , and therefore £ + (^) < rain 2; 



(3.22) 



in the nontrivial case u > k + c, we also have an estimate for the derivative of the logarithmic function: 



\(Q + Y(u)\ 



2-p 



-1 


2-p 


1 


2-p 




< 






+ c + k — u 




C 





( 2 2+n) 



p-2 



(3.23) 



With these estimates at hand, we have for the diffusive term: 

rt 



G 2 ^G\ :=2 (J aiuf-^Vu^Vu- V {g + (u)(g + Y (u)( p } dx ds 

J —i J B R / 2 
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f f a(u) p - 1 \Vu\ p {2(l + g + (u))[(g + )'(u)] 2 C p }dxds + G* 2 as h 0, 

J — t J B R / 2 



where we define 



G\ :=2p [ [ aiuy-^Vuf-^u- VC{^ + (^)(^ + ) / (^)C P " 1 } ^ ^- 

J — i J B R / 2 

Applying Young's inequality (3.8) with the choices 

r=p, a=\Vur 1 C- 1 \(g + y(u)\ 2/p ', b = |(e + ) , («)| 1 " 2/p ' |VC| and e 4 = 1, 

we obtain 

\G* 2 \<2p f I a(M) p - 1 |Vu| p - 1 |VC|^ + (w)|(^ + )'(w)|C p " 1 ^^ 
J-iJB R/2 

= 2pf { a( U f-V + (w)|V U r 1 C p - 1 |(^ + )'(«)| 2/p '|(^ + )'(«)r 2/p '|VCl^^ 

J-tJB R/2 

<2e p f f o(w) J, -V + (w)|Vw| J, [(p + ) , (w)] 2 C p rfa;ds 
+ f f a(ur- 1 g + (u)\VC\ p \(g + Y(u)\ 2 ~ P dxds 

V e 4 V-t JB R/2 

= 2 f f a(w) p -V + («)|Vw| p [(e + )'(u)] 2 C P da;ds 
+ 2(p-l) T / a«-V + («)|VCr|(^ + )'(u)| 2 " P ^d S . 

./— i J B R /2 

In face of this estimate, we obtain 

=2 / a(^- 1 |V^| p [(^ + ) / H] 2 C P dxd5 

J— £ J B R /2 

-2(p-l) /* / a«-V + (w)|VC| p |(^ + )'(«)| 2 " P rfxrf S 

./ — £ J B R /2 

> 2[ 7l ^( W /4)f- 1 f / |v«r[( e + )'w] 2 c p ^^ 

J — £ J Bj} io 



2(p-l) / ^ / a(*if-V(^)|VCr|(£ + )>)| 2 



£ ./£ 



R/2 



and, finally, 



> 2[ 7l ^( W /4)f- 1 /"* /" |V W r[(p + )'(u)] 2 C P ^d S 

J— £ J B R /2 

-2{p- l)nln2 (^T^ / f X ^^X^ 1 ""/ 4 } d;Eds > 

G$ > 2 [ 7l ^ ( W /4)f- 1 /* / |Vw| p [(<?+)'(«)] V <*r ds 

./ — £ J B R /2 

-2(p- l)nln2 (^-) P ~ 2 i|B K/2 | [^(^f- 1 , 



(3.24) 
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where we have used estimates (3.22), (3.23), the properties of £, and the fact that 
7i'0 (cj/4) < a(u) < 72^ (^/4) on the set {u > 1 — cj/4}. 
Moreover, from the definition of t and our choice of t* (recall that t* > dR p — ciqR p ), there holds 

a«RV=("-t P RP r . (3.25) 



t < 



Taking into account (3.25), we obtain from (3.24) that 



Gl > 2 [ 7l V (^/4)] p - 1 f [ I V«|" [(e + )»] V d* 



- 2(p - l)n In 2C2< 1 +-)( a -rt |5* /2 | [72 ^ . 

On the other hand, for the lower order term, by passing to the limit h — > 0, we have 
G 3 ^G*:=2 /"* / xw/(w)VvV W {(l + ^ + (u))[(p + )'( W )] 2 C J '}d a; ds 

^ — £ J B R / 2 

+ 2p ( f X uf(u)Vv-VC{g + (u)(g + Y(u)C p - 1 }dxds 

J —t J Bt> /o. 



(3.26) 



<2M / / (l + ^ + (^))[(^ + ) / (^)] 2 C^|V^||Vv|^ds 



t JB h 



+ 2 P M f [ Q+^Kg+nu)] 1 2/p \vv\\VC\\(Q + )\u)\ 2/p ' C' 1 dxds. 
J-iJB R/2 



R/2 

Applying Young's inequality (3.8) to the first term on the right-hand side with 

r = P , a=\Vu\, b=\Vv\ and e 5 = - - - - j , 

and to the second term with 

r =p, a = |(^+)'(w)| 1 - 2/p ', 6=|Vt;||(p + )'( W )| 2/p 'c p - 1 and e 6 = 1, 

we obtain 

G^<2^( w /4f- 1 



/*/ \Vu\ p [(Q + y(u)] 2 ( p dxds + 2M ff Q+(u)\VC\[(g + Y(u)] 2 p dxds 

J —t J B R / 2 J — t J B R / 2 

+ 2M£ r (stttS)"""' £ + s+w > [<« + )'(«)] Viv.p' 

+ 2M(p-l) / / ^+(u)|VC||Vt;| p '[(^ + )'(zi)] 2 C p rfa;ds. 

J — t J B R / 2 

Using the estimates (3.22) and (3.23) and the properties of f, we then get 

Gl < 2^/4)»" 1 J* J \Vu\P[(Q + nu)] 2 edxds + 2Mnln2^ (^J^ t\B m \ 

+ 2M(p-l)nln2^ (-^) ' / / |V^| p 'x { „>i- w /4 } dxds. 

./ — £ J B R /2 
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Then, applying Holder's inequality and recalling the definition of t, we get 



G% < 2^(uu/Ay- 1 f f \Vu\ p [(Q + )\u)] 2 C P dxds + 2MCn\n22^^ 

J —t J Bf? /o 



+ 2M(p- 1) 



{( 



pip (u/Af 



p-1 x l/(l-p) 



l + nln2 C 



^M(l + nln2)y p 

x \wv\r; p , p(QT) (ao^lB^I) 1 - 1 ^. 

In addition, thanks to Remark 3.1, we may estimate 

V p V(^/4) 



nln2 



2 2+i 



0' 



( 2 2+n) 



M(l + nln2) 



i-i/p 



RV- 1 < 1, 



c(^)"«:-'"«'- 2 <i, *(£)'"■>-' <i, 

and this finally gives 

< 2V'(w/4f- 1 f / |Vw| p [(p + )'(u)] 2 C J 'da; + 2MCn In 2 2( 1+ "» 2 -f)|B fi/2 1 

^ — £ J B R / 2 

+ 2M(p-l)Cnln2||V W ||^, p(Qr) | J B JJ/2 | 1 - 1 ^. 
Combining estimates (3.20), (3.26) and (3.27) yields 

/ {^ + (u)} 2 C p dxd S < 2M(p- \)Cn\n2\\Vvf Lvf \B R/2 \ 



(3.27) 



i-i/p 



(1-7^)2^^/4)] 



p-i 



\\7u\ p [(Q+)'(u)} 2 e dxds 



R/2 



2n In 2 2( 1+ ")( 2 "P) |fl H/2 | J MC + (p - 1)G 



^(w/4) 



(w/2 m )_ 



and since 71 > 1 and n > 0, this implies 

sup / dx < 2M(p- l)Cnln2||V^||^ /pfo .l^l 1 -! 

-t<t<0JB R/2 x{t} VWTJ 



+ 2n\n2 2 2 - p 



\B R /2 I { 



MC + (p- \)C v -il 



p-i 



^)(cj/2 m ) 



p-i 



} 



(3.28) 



Since the integrand in the left-hand side of (3.28) is nonnegative, the integral can be estimated from below 
by integrating over the smaller set S = {x G B R / 2 : u(x,t) > 1 — c<;/2 2+n } C B R / 2 - Thus, noticing that 

C = 1 and {g + (u)} 2 > (ln(2 n - 1 )) 2 = (n - l) 2 (ln2) 2 on S, 
we obtain that (3.28) reads 

\{xeB R/2 : ^(x,t)>l-^/2 2 +-}| 



< 



2Cn\B R/A \ 
(n- l) 2 ln2 



| 2 2-p 



Mc+(p-i)c p 7 r 1 





p-i- 


>(cj/2 m )_ 





M(p-l)||V< p /p 



for all t G (— £, 0). To prove the lemma we just need to choose si depending on v\ such that s\ = 2 + n with 
2C 



n > 1 + 



v x In 2 



MC+(p- l)C p 7 f 



P^p-i 



r " 


p-i- 


_(j)(uu/2 m )_ 





M(p-l)\\W\\l P , p 



Lp'p(Qt) [' 
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since if n > 1 + 2/a then n/(n — l) 2 < a, a > 0. Furthermore, si is independent of uo because 



r j 


p-i - 


>(^/2 m )_ 





(w/4) 



/3 2 /b-i) 



(u;/2 m ) 



/3i/b-i) 



(p-i) 



,/3 2 -/3i 9 m/3 1 -2/3 2 



< 2 



m/3!-2/3 2 



The last inequality holds since fa > Pi- 

Now, the first alternative is established by the following proposition. 



□ 



Proposition 3.1. The numbers v\ G (0,1) and s\ ^> 1 can 6e chosen a priori independently of uo and R, 
such that if (3.4) holds, then 



2 si+i 



a.e. in Q(t,R/8). 



We omit the proof of Proposition 3.1 because it is based on the argument of [ , Lemma 3.3] and [7], and 
we may use for the extension the same technique applied in the proof of Lemma 3.2. 

Corollary 3.1. There exist numbers vq,o~q G (0, 1) independent of uo and R such that if (3.4) holds, then 

essosc u < cfquo. 

Q(i,R/8) 

Proof In light of Proposition 3.1, we know that there exists a number si such that 



ess sup u < 1 

Q(i,R/8) 



2 fi i+ i: 



and this yields 



ess osc u = ess sup u — ess inf u < I 1 



2*1+1 



Q(t,R/8) Q(f,R/8) Q^R/S) 

In this way, choosing ctq = 1 — l/2 Sl+1 , which is independent of we complete the proof. 



□ 



3.4. The second alternative. Let us suppose now that (3.4) does not hold. Then the complementary case 
is valid and for every cylinder we have 



\{(x,t) G : u(x,t) <uj/2}\ < (l-i/ )|Qg|. 



(3.29) 



Following an analogous analysis to the performed in the case in which the solution is near its degeneracy 
at one, a similar conclusion is obtained for the second alternative (cf. [4] and [ ]). Specifically, we first 
use logarithmic estimates to extend the result to a full cylinder and then we conclude that the solution is 
essentially away from in a cylinder Q(r, p). In this way we prove the following corollary. 

Corollary 3.2. Lett denote the second- alternative- counterpart oft. Then there exists g\ G (0, 1) ; depending 
only on the data, such that 

essosc u < <j\u. 

Q(t,R/8) 

Since (3.4) or (3.29) must be valid, the conclusion of Corollary 3.1 or 3.2 must hold. Thus, choosing 
o~ = max{<To,<7i} and t° = min{t,t}, we obtain the following proposition. 

Proposition 3.2. There exists a constant a G (0, 1), depending only on the data, such that 

essosc u < auj. 

Q(t<>,R/8) 



The local Holder continuity of u in Q T now follows (see, e.g., [ ], [ ], or the proof of [23, Th. 2]). 
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Figure 1. Example 1: Numerical solution for species u, at t = 1.0 for p = 2 (left), and p = 6 (right). 



4. Numerical examples 

In this section, we provide two numerical examples to illustrate how the approximate solutions of the 
chemotaxis model (1.1) vary when changing the parameter p from standard nonlinear diffusion (p = 2) to 
doubly nonlinear diffusion (p > 2). For the discretization of both examples, a standard first order finite 
volume method (see the Appendix for details on the numerical scheme) on a regular mesh of 262144 control 
volumes is used. We choose a simple square domain Q = [— 1, l] 2 and use the functions a(u) = eu(l — u), 
f(u) = (1 — u) 2 and g(u, v) = au — (3v, along with parameters that are indicated separately for each case. 



4.1. Example 1. For the first example, we choose e 
initial condition for the species density is given by 



0.01, a = 40, (3 = 160, x 0.2 and d = 0.05. The 



1 for \\x\\ < 0.2, 
otherwise, 



and the chemoattractant is assumed to have the uniform concentration ^o(^) =4.5. In a first simulation, we 
consider the simple case of p = 2 and we compare the result with an analogous experiment with p = 6. We 
evolve the system until t = 1.0, and show in Figure 1 a snapshot of the cell density at this instant for both 
cases. 



4.2. Example 2. We now choose the parameters e 
condition for the species density is given by 



u Q (x) = 



0.5, a = 5, f3 = 0.5, x = 1 and d = 0.25. The initial 

1 for \\x- (-0.25,0.25)|| < 0.2 or \\x - (0.25, -0.25)|| < 0.2 
otherwise, 
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Figure 2. Example 2: Numerical solution for species u, at £ = 0.1 for p — 2 (left), and p = 6 (right). 

and for the chemoat tract ant 

v ( x \ _ j 45 for (0.25,0.25)|| < 0.2 or \\x + (0.25, 0.25)|| < 0.2 
1 otherwise. 

The behavior of the system for the cases p = 2 and p = 6 at different times is presented in Figures 2, 3 and 4. 

4.3. Concluding remarks. W first mention that, from the previous examples, one observes that even 
though the numerical solutions obtained with p = 2 differ from those obtained with p > 2, the qualitative 
structure of the solutions remains unchanged. We also stress that the numerical examples illustrate the 
effectiveness of the mechanism of prevention of overcrowding, or volume filling effect, since all solutions 
assume values between zero and one only. In particular, all examples exhibit plateau-like structures where 
u = u m = 1, at least for small times, which diffuse very slowly, illustrating that the diffusion coefficient 
vanishes at u = 1 (recall the special form of the functions a(u) and f(u): they include the factor (1 — u), 
and therefore the species diffusion and chemotactical cross diffusion terms vanish at u = and u = u m = 1). 

In Example 2, the solution for p — 2 has a smoother shape than the one for p = 6, which exhibits sharp 
edges. These sharp edges do not only appear for u = and u = u m , where one expects them, due to the 
degeneracy of the diffusion term and the choice of initial data, but also for intermediate solution values, as 
is illustrated by the plots for p = 6 of Figures 2 and 3. 
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Figure 3. Example 2: Numerical solution for species u, at t = 0.5 for p = 2 (left), and p = 6 (right). 
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Appendix 



The definition of the finite volume method is based on the framework of [28]. An admissible mesh for Q 
is given by a family T of control volumes of maximum diameter h, a family of edges £ and a family of points 
(xx)KeT- For K G T, xk is the center of K, £[ nt (K) is the set of edges a of if in the interior of T, and 
£ ext (K) the set of edges of K on the boundary dCl. For all a G £, the transmissibility coefficient is 



d(x K ,x L ) 

M 



for a e £int(K), cr 
for a e £ext(K), 



K\L, 



K d(x K ,a) 

where K\L denotes the common edge of neighboring finite volumes K and L. For K E T and cr = K\L £ £{K) 
with common vertexes (cbi,K,L)i<i<i with / G N\{0}, let T G (Tj^ for cr G £ ext (if ), respectively) be the open 
and convex polygon built by the convex envelope with vertices (xk,xl) {xr, respectively) and (cli^k,l)i<£<i- 
The domain Q can be decomposed into 
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Figure 4. Example 2: Numerical solution for species u, at t = 2.5 for p = 2 (left), and p = 6 (right). 



For all K G T, the approximation X7hUK,a of Vii is defined by 

if a = K|LGf int (X), 







ifdGfextW- 



To discretize (1.1), we choose an admissible mesh of Q and a time step size At > 0. If Mt > is the smallest 
integer such that M T At > T, then t n := nAt for n G {0, . . . , M T }. 

We define cell averages of the unknowns A(u), f(u) and g(u,v) over K e T : 

A(ia(x, £)) dt, ^ 7 ^ rl := ~^t\K\ I I v ( x ^)) d% dt, 



An+l — x 

Ar - m\k\ 



K 



K 



fn+1 ._ 
JK 



1 



At I if I 



K 



f(u(x, t)j dx dt, 



and the initial conditions are discretized by 

U ° K = W\ Ik V ° K = W\ Ik 

We now give the finite volume scheme employed to advance the numerical solution from t n to t n+1 , which is 
based on a simple explicit Euler time discretization. Assuming that at t = t n , the pairs (u K , v 7 ^) are known 
for all K G T, we compute (u 7 ^ 1 , v 7 ^ 1 ) from 



7y n+l 

\k\ Uk 



l K 



At 



E r a \V h A K X 2 ^* A 



aeS(K) 



aeS(K) 



(V h v^) + u K fl-{V h v^ 
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\K\ K M K = £ T^ h v^ a + \K\gl. 

Here | • |^ denotes the discrete Euclidean norm. The Neumann boundary conditions are taken into account 
by imposing zero fluxes on the external edges. 
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